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Abstract 

The many experimental puzzles of solid hydrogen, that seem to defy the generally 
accepted view of condensed matter, are shown to be a rather straightforward con- 
sequence of a new approach to condensed matter physics, that takes into account 
the coherent electrodynamic processes, predicted by Quantum ElectroDynamics 
(QED), that have been so far overlooked. As a consequence, a revision of the 
present view, based on electrostatic forces only, is called for. 
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Much as the hydrogen atom has been the paradigm-problem of Atomic Physics, 
it would seem that solid hydrogen should have the same role in solid state physics, 
a simple and essential testing ground of the prevailing theoretical ideas on con- 
densed matter. And yet, solid hydrogen has become over the years an embarassing 
puzzle, a source of remarkable discoveries that have found the generally accepted 
paradigrruthat somewhat impressionistically we may describe as the Electrostatic 
Meccano lil, quite off-balance. Indeed, there are two major aspects of solid hydro- 
gen phenomenology that, in our opinion, expose the inadequacy of the present 
view: at low pressure the observed equation of state defies all efforts based on 
"reasonable" two-body potentials, and appears to require the action of "many- 
body forces" 0, while at high pressure metallic hydrogen, predicted by the classic 
work of Wigner and Huntington || to occur at pressures of 25 GPa, at pressures 
ten times larger still eludes detection. On the other hand, for P ~ 150 GPa 
a new remarkable phase of molecular hydrogen, namely solid hydrogen III, has 
been discovered about ten years ago, characterized by intense infrared absorption 
and a sharp discontinuity in the vibron frequencies H, [ffl. 

In this Letter we wish to report a few results on the physics of solid molecular 
hydrogen of a new approach to condensed matter which, based on QED, takes 
finally into account the non-negligible collective electrodynamic actions that the 
elementary systems (atoms and molecules) exert upon each other at high enough 
densities and low enough temperatures. As the basic theory, as well as a few 
relevant applications to the physics of condensed matter have been published m 
a recent, readily accessible book ||, here we shall refrain to review the theory @, 
and only mention the two theoretical developments which, as we will show, have 
an important bearing upon the physics of solid molecular hydrogen, namely the 
application of the QED coherent interaction to the solidification of 4 He [|7| 
and to the dynamics and thermodynamics of liquid water ||, H. 

Let us then address first the problem of solidification of molecular hydrogen 
at low pressures. We shall follow London |J| in modeling the solid as an ensemble 
of cages where each molecule "rattles", confined by its nearest neighbours. As we 
are not interested in a detailed calculation of the spatial structure of the solid, but 
only on the collective electrodynamic interactions, that have been so far neglected, 
we shall model the cage as a spherical cavity [0 of radius R = r — r*, where r 
is the average intermolecular distance and r* is the radius of the co- volume, that 
we evaluate by imposing that at r* the total energy, kinetic plus Lennard- Jones, 
be zero. For H 2 and D 2 we obtain r* = 2.79 A and 2.80 A respectively. As 
explained in Refs.0, 0, the main idea is that due to its repulsive interaction 
with the nearest neighbours, the walls of its cage, each molecule, slightly violating 
the Born-Oppenheimer approximation, develops a weak electric dipole, through 
which it gets coupled to the quantized modes of the electromagnetic field. By 
means of these simple considerations we may look at our molecules as two-level 
systems, "oscillating" between the s-wave ground state and the p-wave first- 
excited state, whose energies are 

f - 1x2 e - (449)2 m 



1 According to the generally accepted view, the interactions among the atomic/molecular 
systems of condensed matter have an electrostatic nature, hence a short-range character, just 
like the mechanical forces that hold together the pieces of a Meccano. 

2 We believe it should also be of some interest to the reader to consult a recent article by 
CP. Enz, where the basic theoretical ideas of QED coherence are critically analysed. 
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respectively I, resonantly coupled to the electromagnetic modes of frequency 

u-E E - 10 ' 29 (2) 
2m(r — r*) 2 

In order to complete the quantitative characterization of the two-level system we 
must estimate the size of its electric dipole, which determines the strength of its 
coupling to the radiation field. As argued in Refs. ||, 0, a reliable estimate 
of the dipole D{r) = 2ed{r) can be obtained by equating the repulsive potential 
of two dipoles at the distance r with the repulsive part of the Lennard- Jones 
potential, i.e. by setting 

Af^r^, ,3) 



where A=36.7 K and i? m =3.32 A for H2 (35.2 K and 3.31 A in D2) as determined 
by the second virial coefficient of the gas phase [|l(|, and a = is the fine 
structure constant. 

Once we know u>(r), and the effective coupling constant given by 

. (Ak\ 1 / 2 2ed(r) (N\ 1 ' 2 
9{r)=[ — ) -7±k 77 , (4) 
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the coherent many-body electrodynamic interaction is fully determined f|, 0], 
and we find that for P = and T = 0, at the density of the solid 

/16\ 1/2 

g (r = 3.79A) = 1.72 > g c = (-J ^0.77, (5) 

Similarly for deuterium we find g(r = 3.7lA)=2.32, that is we are in the regime 
of strong coherence for both isotopes. In Figure 1 we report the energetics of 
both solid hydrogen and deuterium as a function of the intermolecular distance 
r. Please note that uo{r) — > 00 as r — > r* and Eq.(4) predicts that g{r) < gc, 
thus for such values of r the energy gain from the long-range coherent interaction 
disappears. For the zero-pressure binding energies of a molecule in the lattice, 



we obtain the values E— - 83 K for H 2 (exp. - 85.5 K |TT|) and E— - 127 K for 
D 2 (exp. -132 K [jH|). It should also be perfectly clear that for this mechanism to 
be operative one needs a well defined u>(r), a requirement that only a crystalline 
structure can meet. 

Without much ado, and in an obviously perfectible way, we may extend our 
calculation for T 7^0. In Ref. it is shown that the disordering effect of thermal 
fluctuations produces a two-fluid system, comprising an ordered solid in patches 
of the size of the Coherence Domain (CD), of radius Rqd — ^y, and a disordered 

fluid that resides in the interstices. The T-dependence of the coherent, ordered 
fraction is given by (Figure 2): 



Pc {T) = 1 - J* C ° r 2 exp [--^) dr, (6) 



3 We are using here the "natural units" , where h = c — fcs = 1 
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where Vqd is the volume of the CD, and 5(r) is a well defined r-dependent gap ||, 
that arises from the coherent electrodynamic process we have described above. 
Thus, for the free energy per molecule in the solid we may write: 



E e s Hq U id -\- 5, 



F solid {T)=p c {T)8 + {l-pc{T)) 

L \ \ X V / so lid / J 

( 7 ) 

where Z rot is the rotational partition function of the molecule, and 8 e _ s _ is the 
"electrostatic" gap between the solid and the liquid phase, which in the mean 
field approximation, by use of the pair-potential and the values of the molar 
volumes in the liquid Vl = 26.1 cm 3 and in the solid Vs = 23.15 cm 3 , turns out 
to equal 5 e , s . = E e . s ., S oiid — E e , s .,iiq U id — —4 K. In a similar fashion we have for the 
liquid free energy per molecule: 



liquid 



(T) = -Tin 
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E es . : liquid (8) 



where Z trans is the translational partition function. In Figure 3 we have plotted 
the free energy difference per molecule between the solid and the liquid phase, 
as a function of temperature for P—0. We see that the two curves cross zero at 
T=13.5 K for H2 and T—17 K for D2, in adequate agreement with the observed 
solidification temperatures Tsx=13.96 K of H 2 and Ts£=18 K of D 2 respectively. 

For finite T, as the pressure increases, London's spherical cage ceases to be 
a good approximation, as we can easily understand from the vanishing of the 
gap when r — > r* (see Fig. 1). Thus at high pressure we expect the precise 
spatial structure and the rotational state of the molecules to realize a new " cage 
configuration", whose energetic advantage will possibly lead, as observed, to a 
transition to the phase II. 

Leaving the complex and fascinating subject of phase II at this qualitative 
level of understanding, we wish now to address the extremely intriguing problem 
of phase III, reached at pressures of about 150 GPa. As well known, the signal 
of this phase transition is a gigantic increase in the absorbance of the infrared 
(IR) vibron mode and a sharp discontinuity Av in the frequencies of both IR and 
Raman vibrons. Both facts have been, in our opinion, correctly interpreted in 
a recent paper JES, as arising from a sudden change in the electronic configura- 
tion of the H2 (JJ2) molecule, which thus acquires a permanent electric dipole. 
However, we think the mechanism of dielectric instability, proposed in the above 
work to account for this most interesting phase transition, hardly believable. We 
shall now show that QED coherence has something very relevant to say also for 
phase III. Indeed, even though the electromagnetic couplings of the molecular 
ground state to the excited states (AE > 10 eV) is known to be rather small, 
like it happens in the H2O molecule f|, the access to a coherent ground state is 
made possible by a photon "mass term" inside the solid given by: □ 



2 uj(E n -E a ) 

V V / n \tj n — hi a ) - U- ij =ll2 ,3 r=l,2 



where k refers to the electromagnetic field mode with momentum k, uj^ is the 

energy of the mode (oj? = \k\ in free space) and \n) is a complete set of molec- 
ular states of energy E n . By performing an analysis totally similar (but for the 



4 See for instance Ref. H, pages 50-53 and 195-205 
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different experimental inputs: oscillator strengths and energy spectrum) to the 
one carried out for water [pfl, M, we obtain the following results: 



the threshold for the coherent electrodynamic process that rearranges 
the electronic configuration of the H2 {D2) molecule occurs at 



P_ 

Po 



(10) 



to be compared with ( — ) ~ 9 

V PO / exp 

in the new configuration the molecules are in a state 

\C) = cose\0) +sm6e- luJrt \U u ), (11) 

with 

sin = 0.28, (12) 

a superposition of the ground state |0) and the molecular configuration 
whose energy in the gaseous state is E(H U )—12A eV; 

the energy gain per molecule, the gap S e , turns out to be rather small: 

5 e ~ - 0.10 eV. (13) 
Estimating the electric dipole of the \U U ) state by analogy with the HC1 



molecule [13 



— ~ 2 ■ 10~ 9 cm (14) 

e 

yields an oscillator strength for the IR vibron 

2 / D\ 2 

fth ^ -m e u vibron sin 4 9 f — J = 1.13 • 10" 5 , (15) 

to be compared with f exp ~ 1.2 • 10~ 5 [Q]. 

With the above informations we may also evaluate the frequency shifts Au 
of the vibrons in Phase III. Indeed, our model predicts that, as the threshold 
for electronic coherence is reached, the vibrational dynamics of the protons will 
undergo a discontinuous change. In fact, when the electronic wave function ac- 
quires a small but finite admixture of the dipolar configuration |II U ), the protons 
no longer oscillate in the previous potential well. Quantitatively, we took the 
experimentally known frequencies of the intramolecular stretching modes in the 
gaseous phase for the ground and first-excited electronic levels , in which the 
vibron has the values i/=4395 (3118) cm" 1 and u*=2U2 (1736) cm" 1 in H 2 (D 2 ) 
respectively | |T3|] , and the value of the "mixing angle" 9 ~0.29. 



In the coherent electronic ground state (Eq. 11) which marks phase III, the 
old vibron is replaced by a new vibron, whose frequency is: 

via — v cos 2 + v* sin 2 9 (16) 
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i.e. 4235 cm -1 in H2 and 3005 cm -1 , yielding downward shifts Au = 160 cm -1 in 
H2 and 113 cm -1 in D2, only somewhat larger than those observed by Mao and 
Hemley B at finite (T=77 K) temperature, i.e. about 103 cm -1 for deuterium. 
In view of our neglect of the effect of pressure and temperature on the vibron 
frequencies, we regard this result as quite satisfactory. 

In Figures 4a, 4b we present our results for the equation of state of solid hydro- 
gen compared with recent experimental data |]T3J. It provides better agreement 
than other calculations based on a purely electrostatic approach [16] . 



In conclusion, we believe we have provided sufficient evidence to recommend a 
serious and rational reconsideration of the prevailing views of condensed matter, 
especially in the light that the approach, which we have shown to account effort- 
lessly for many puzzles of solid hydrogen, encompasses all the QED interaction 
mechanisms, not merely the electrostatic ones. 
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